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1.  INTRODUCTION:  This  project  involves  the  development  of  optically  guided  multi-agent 
therapies  for  metastatic  and  drug  resistant  breast  cancer.  Non-invasive  Diffuse  Optical  Imaging 
technologies  are  able  to  monitor  drug  response  and  resistance  through  quantitative  tracking  of 
tumor  metabolism  and  vascular  supply,  using  clinic  friendly  and  portable  devices  that  can 
monitor  deep  lesions  (multiple  centimeters  below  the  skin)  in  both  breast  tissue  and  at  the  site  of 
bone  metastases.  To  demonstrate  our  strategy,  we  will  show,  for  the  first  time,  that  non- 
invasive  optical  guidance  can  enhance  the  efficacy  of  combined  cytotoxic  and  antiangiogenic 
therapy  while  simultaneously  increasing  the  time  until  relapse  and  decreasing  toxicities 
associated  with  high  drug  dosing.  This  will  be  accomplished  without  the  use  of  exogenous 
contrast  agents  or  dedicated  imaging  facilities.  Additionally,  in  order  to  demonstrate  the 
feasibility  of  translating  optical  signatures  discovered  in  animal  models  to  breast  cancer  patients, 
we  will  conduct  a  clinical  study  of  non-invasive  optical  monitoring  of  drug  response. 

2.  KEYWORDS:  drug  resistance,  imaging,  optics,  diffuse  optical  imaging,  spatial  frequency 
domain  imaging,  intravital  microscopy,  multiphoton  imaging,  therapy  monitoring,  metastases, 

3.  ACCOMPLISHMENTS: 

■  What  were  the  major  goals  of  the  project? 

■  For  reference,  the  original  SOW  table  is  shown  below,  completed  Tasks  and  Subtasks  are 
shaded  in  dark  blue,  ongoing  tasks  are  in  lighter  blue. 


Specific  Aim  1:  Preclinical  Investigation  of 

Systemic  Therapies  with  Optical  Imaging 

Timeline 

Site  1 

Major  Task  1:  Instrument  Setup  and  Testing 

Months 

Subtask  1 :  Local  IACUC  Approval 

1-3 

Dr.  Roblyer 

Subtask  2:  Milestone  Achieved:  ACURO 
Approval 

3-6 

Dr.  Roblyer 

Milestone  #1  ACURO  approval  obtained 

6 

Dr.  Roblyer 

Subtask  3:  Custom  SFDI  Fabrication,  order 
parts  and  construct  device 

1-9 

Dr.  Roblyer 

Subtask  4:  Setup  Multiphoton  Microscope, 
work  with  vendors  to  order  and  install/test 
microscope 

1-9 

Dr.  Roblyer 

Subtask  5:  Perform  initial  SFDI  and  MPM 
animal  testing  with  SHO  mice  and  MDA-MB- 
231  cells. 

[12  mice  x  1  groups  =  12  mice] 

9-12 

Dr.  Roblyer 

Subtask  6:  Evaluate  growth  rates  and  treatment 
response  for  SHO  mice  and  MDA-MB-231  cells. 
[12  mice  x  1  groups  =  12  mice] 

6-12 

Dr.  Roblyer 

Milestone  #2  SFDI  and  MPM  ready  for  animal 
imaging 

12 

Dr.  Roblyer 

Major  Task  2:  Cytotoxic  +  Antiangiogenic 
Monitoring 

Subtask  7:  Evaluate  treatment  timelines  and 
response  characteristics  for  SHO  mice  and 
MDA-MB-231  cells/  Monitor  with  SFDI. 

[18  mice  x  3  groups  =  54  mice] 

12-24 

Dr.  Roblyer 

Subtask  8:  Correlate  intravital  MPM  with  SFDI 
[10  mice  x  1  groups  =  10  mice] 

12-24 

Dr.  Roblyer 

Subtask  9:  Test  growth  rates  of  MMTV-PyMT 
mice 

[10  mice  x  1  groups  =  10  mice] 

12-24 

Dr.  Roblyer 

Milestone  #3  Define  optical  signatures  of 
cytotoxic  and  antiangionic  therapy  response 

24 

Dr.  Roblyer 

Milestone  #4  Define  quantitative  correlates 
between  MPM  vascular  imaging  and  SFDI 

24 

Dr.  Roblyer 

Major  Task  3:  Optically  Defined  Cytotoxic  + 
Antiangiogenic  Therapy 

Subtask  10:  Test  scheduling  of  cytotoxic  + 
antiangiogenic  based  on  optical  signatures, 
(dependent  on  Major  Task  2) 

[18  mice  x  3  groups  =  54  mice] 

24-36 

Dr.  Roblyer 

Subtask  1 1 :  Additional  correlation  of  intravital 
MPM  and  SFDI  in  optical  defined  therapeutic 
window. 

[10  mice  x  1  groups  =  10  mice] 

24-36 

Dr.  Roblyer 

Subtask  12:  Test  SFDI  and  MPM  imaging  with 

MMTV-PyMT  mice 

[10  mice  x  1  groups  =  10  mice] 

24-36 

Dr.  Roblyer 

Milestone  #5  Improved  survival/outcomes 
though  optically  defined  therapy  scheduling 

36 

Dr.  Roblyer 

Major  Task  4:  Early  Response  Monitoring 

Subtask  13:  Test  early  (hrs)  response  with  SFDI, 
MPM,  IHC,  determine  immune  modulators 
(dependent  on  Major  Task  2  and  3  results) 

[18  mice  x  3  groups  =  54  mice] 

36-48 

Dr.  Roblyer 

Milestone  #5  Characterize  immune  modulators 
of  early  optical  response  to  systemic  therapy. 

48 

Dr.  Roblyer 

Major  Task  5:  New  Therapy  Monitoring 

Subtask  14:  Test  optical  signatures  of  response  to 
her2/hormone  and/or  immunotherapies. 

(dependent  on  Major  Task  2,3  and  4  results) 

[18  mice  x  3  groups  =  54  mice] 

48-60 

Dr.  Roblyer 

Milestone  #6  Define  optical  response  of  new 
therapies. 

60 

Dr.  Roblyer 

Specific  Aim  2:  In-vivo  Clinical  Studv  of 

Progressive  Resistance 

Major  Task  6:  dDOS  fabrication 

Subtask  15:  Design/Fabricate  dDOS  system  and 
new  custom  dDOS  probe 

6-24 

Dr.  Roblyer 

Milestone  #7  Clinical  Ready  dDOS  system 

24 

Dr.  Roblyer 

Major  Task  7:  Normal  Volunteer  Study 

Subtask  15:  Local  IRB  Approval  (for  both 
normal  volunteers  and  clinical  study) 

30-33 

Dr.  Roblyer 

Subtask  16:  HRPO  Approval 

33-36 

Dr.  Roblyer 

Subtask  17:  Measure  10  normal  volunteers  with 
appropriate  updates  to  probe.  Analyze  and 
interpret  imaging  data. 

[10  normal  volunteers] 

36-48 

Dr.  Roblyer 

Milestone  # 8  Normal  Volunteer  Study 

Completed  and  dDOS  system/probe  ready  for 
clinical  study. 

48 

Dr.  Roblyer 

Major  Task  4:Breast  Cancer  Clinical  Study 

Subtask  18:  Measure  30  breast  cancer  patients. 

[30  breast  cancer  patients] 

36-57 

Dr.  Roblyer 

Subtask  19:  Analyze  and  interpret  imaging  data. 

36-60 

Dr.  Roblyer 

Milestone  #9  Clinical  Study  Completed,  Results 
Published 

60 

Dr.  Roblyer 

■  Completion  dates  for  subtasks  are  listed  here: 


■  Subtask  1 :  Local  IACUC  Approval:  5/31/2016 

■  Subtask  2:  ACURO  Approval:  8/1 7/201 6 

■  Subtask  3:  Custom  SFDI  Fabrication:  ongoing 

■  Subtask  4:  Setup  Multiphoton  Microscope:  8/1/201 6 

■  Subtask  5:  Initial  SFDI  and  MPM  testing:  7/1/2017 

■  Subtask  6:  Evaluate  growth  rates:  7/1/2017 

■  Subtask  7:  Treatment  monitoring  with  SFDI:  ongoing 

■  Subtask  8:  Correlate  intravital  MPM  with  SFDI:  ongoing 

■  Subtask  9:  Evaluate  growth  rates:  7/1/2017 


■  What  was  accomplished  under  these  goals? 

■  Major  accomplishments  and  results  are  listed  for  each  project  goal  for  year  2  as  stated 
in  the  subtasks  in  the  SOW.  Note:  updates  from  year  1  for  completed  tasks  are  not 
included  in  this  report. 


■  Subtask  3:  In  prior  reports,  updates  on  custom  SFDI  Fabrication,  including  new  angle 

correction  algorithms  and  other  technical  updates  were  provided.  Here  we  provide  additional 
updates.  It  is  of  note  that  this  subtask  has  become  a  fruitful  instrumentation  development 
aspect  of  the  project,  and  will  likely  continue  throughout  the  award  period. 

■  A  manuscript  describing  a  new  way  to  conduct  simultaneous  multi-wavelength 
(hyperspectral)  Spatial  Frequency  Domain  Imaging  (SFDI)  has  been  published 
(see  appendix). 

■  SFDI  short-wave  infrared  instrument  development:  SFDI  has  been  developed 
to  measure  wavelengths  spanning  from  the  near-infrared  to  the  short-wave 
infrared  (approx.  690nm  -  1300nm).  This  allows  for  the  extraction  of  quantitative 
chromophore  concentrations  of  oxy  and  deoxyhemoglobin,  water,  and  lipids. 
These  parameters  are  highly  relevant  to  oncology.  A  manuscript  is  current  in 
progress  describing  this  new  technique.  Below  is  a  diagram  of  the  system  along 
with  imaging  processing  steps. 


Absorption  Scattering 


The  figure  below  shows  SWIR  SFDI  data  collected  from  resected  mouse  tumors,  showing  that 
water  and  lipids  can  be  mapped  over  the  cut  tumor  surface.  This  is  novel,  and  will  likely  lead  to 
new  ways  to  quantify  tumor  heterogeneity  in  both  the  preclinical  and  clinical  settings. 


SWIR  SFDI  -  Tumor  Imaging 


Subtask  6:  Evaluate  growth  rates  and  treatment  response  for  SHO  mice  and  MDA-MB-231  cells. 

■  Initial  in  vitro  IC50  cytotoxicity  testing  has  been  conducted  for  two  MMTV-PyMT 
breast  cancer  cell  lines.  As  described  in  detail  in  section  5,  we  have  now  decided 
to  use  MMTV-PyMT  derived  cell  lines  rather  than  the  MDA-MB-231  cell  lines. 

We  have  found  that  the  Py230  cell  line  (triple  negative,  murine,  mesenchymal) 
grow  too  rapidly  in  vivo  and  form  scars  and  have  a  very  poor  take-rate 
(approximately  50%).  We  found  that  the  Py81 19  cell  line  (murine,  epithelial),  are 
slow-growing,  do  not  form  scars,  and  have  almost  a  100%  take  rate.  The 
decision  was  made  to  move  forward  with  the  use  of  this  cell  line  for  portions  of 
the  future  work. 

Subtask  7:  Evaluate  treatment  timelines  and  response  characteristics  for  SHO  mice  and  MDA- 

MB-231  cells/  Monitor  with  SFDI. 

■  This  subtask  is  ongoing  and  may  involve  the  use  of  Herceptin  responsive  and 
Herceptin  resistance  cell  lines. 

Subtask  8:  Correlate  intravital  MPM  with  SFDI 

■  Correlations  between  microvascular  architecture  (e.g.  vascular  density, 
tortuosity,  etc.)  imaged  using  intravital  MPM  and  global  hemoglobin  levels  are 
currently  being  analyzed  in  a  cohort  of  mice  with  Py230  tumors.  A  figure 
showing  this  type  of  imaging  and  analysis  is  shown  below.  A  manuscript  is 
currently  in  preparation  to  describe  this  multiscale  imaging  procedure. 


Goal  1:  Develop  Multiscale  Widefield  Diffuse 
and  Multiphoton  Imaging  Platform 

Widefield  Spatial  Frequency  Multiphoton 
Domain  Imaging  Microscopy 


Chamber  Model  Tumor:  131  pM 

Subtask  9:  Test  growth  rates  of  MMTV-PyMT  mice 

■  The  growth  rates  Py81 19  cell  line  (murine,  epithelial)  have  been  characterized. 
Subtask  15:  Design/Fabricate  dDOS  system  and  new  custom  dDOS  probe 

■  Circuits  have  been  designed  and  printed  circuit  boards  are  currently  being  laid 
out.  New  laser  sources  for  the  clinical  system  are  being  tested.  New  detector 
modules  are  being  designed  and  fabricated. 


p-vessel 

density 


p-vessel 

tortuosity 


■  Additional  Items  not  explicitly  covered  in  Subtasks: 

■  Significant  training  and  protocol  optimization  has  been  conducted  for 

immunohistochemistry  (IHC).  These  procedures  will  be  important  as  the  project 
goes  forward  as  many  imaging  parameters  will  require  validation  with  IHC  and 
molecular  markers  tissue  testing. 


■  What  opportunities  for  training  and  professional  development  has  the  project  provided? 

■  All  Ph.D.  students  and  postdocs  attend  weekly  lab  meetings  where  they  present  their 
research  results  and  have  regular  meeting  with  the  Ph.D.  Advisor  (Roblyer) 

■  All  Ph.D.  students  and  postdocs  attend  Journal  Club  meetings  where  a  different  student 
presents  a  relevant  paper  in  breast  cancer  and/or  biomedical  optics,  accompanied  by  group 
discussions. 


■  Several  Ph.D.  students  attended  or  gave  presentations  at  research  conferences  over  the  past 
year  as  described  in  section  6. 

■  The  PI  (Roblyer),  gave  multiple  invited  seminars  listed  in  section  6. 

■  How  were  the  results  disseminated  to  communities  of  interest? 

■  Results  were  disseminated  through  peer-reviewed  publications,  conference 
presentations,  and  invited  seminar,  listed  under  part  6:  Products. 

■  What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

The  next  major  push  for  the  project  is  to  utilize  the  new  imaging  technology  development  to 
start  to  define  optical  signatures  of  chemo-response  and  chemo-resistance,  as  well  as  define 
how  different  therapy  schedules  elicit  different  biological  responses  (immunologic  versus 
antiangiogenic).  In  order  to  do  this,  during  the  next  reporting  period,  we  plan  to  complete 
subtask  8,  culminating  in  a  publication  (in  preparation),  continue  work  on  Subtask  3  with 
publication  (in  preparation),  continue  work  or  start  work  on  subtasks  7,10,11,12. 

4.  IMPACT:  Describe  distinctive  contributions,  major  accomplishments,  innovations,  successes,  or  any 
change  in  practice  or  behavior  that  has  come  about  as  a  result  of  the  project  relative  to: 

■  What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

■  The  development  of  SWIR  SFDI  is  likely  to  have  an  impact  on  the  field  of  diffuse  optical  and 
preclinical  oncology  imaging,  as  deeper  tissue  penetration  as  well  as  the  ability  to  quantify 
new  chromophores,  including  lipids  and  collagen,  may  be  highly  relevant  to  chemotherapy 
and  resistance  monitoring. 

■  The  development  of  hyperspectral  multi-wavelength  SFDI  will  likely  impact  both  preclinical 
and  clinical  diffuse  optical  procedures  as  our  recently  published  procedure  is  rapid  and  is  not 
sensitive  to  ambient  lighting  conditions. 

■  What  was  the  impact  on  other  disciplines? 

■  As  the  project  is  still  early,  there  is  nothing  to  report  at  this  time. 

■  What  was  the  impact  on  technology  transfer? 

■  Nothing  to  Report. 

■  What  was  the  impact  on  society  beyond  science  and  technology? 

■  The  project  has  the  aim  of  impacting  society  by  improving  the  health  and  survival  of  breast  cancer 
patients. 


5.  CHANGES/PROBLEMS: 


Changes  in  approach  and  reasons  for  change 

This  change  was  previously  included  in  past  quarterly  reports:  Additional  literature  review  and 
discussions  with  other  researchers  have  taken  place  relation  to  the  choice  of  animal  models. 
In  the  last  progress  report  we  had  decided  to  alter  the  animal  models  from  the  MDA-MB-231 
(cell  line)  SCID  Hairless  Outbred  (SHO)  mouse  model  to  a  syngeneic  mammary  mouse 
model  (Balb/c  mouse)  with  either  4T 1  or  EMT6  mouse  mammary  tumor  cell  lines.  The  cell 
lines  are  current  being  tested  for  growth  rates  and  chemosensitivity  using  MTT  assay. 
However,  we  have  learned  that  both  the  4T1  and  EMT6  are  fast  growing  and  often  develop 
scars  on  the  skin  covering  the  tumor,  which  is  a  major  barrier  for  animal  imaging.  In  addition 
to  these  two  cell  lines  we  have  now  investigated  two  additional  cell  lines  derived  from  the 
MMTV-PyMT  spontaneous  breast  mouse  model  to  determine  if  their  growth  characteristics 
are  more  favorable.  One  of  these  cell  lines,  the  Py230,  has  favorable  characteristics  and 
is/was  used  for  subtasks  5,  6,  7,  8,  and  9.  This  is  a  change  in  the  specific  cell  line,  but  not  in 
the  overall  project  goals. 

For  subtasks  7,10,11,12  and  13,  the  goals  of  these  tasks  has  been  further  elucidated.  For 
subtasks  7,10,11,12  we  will  define  optical  signatures  of  Herceptin  chemo-response  and 
chemo-resistance.  For  subtask  13  we  will  define  how  different  therapy  schedules  elicit  very 
different  biological  responses  (immunologic  versus  antiangiogenic).  These  goals  are 
consistent  with  the  original  proposal  plan. 

Although  not  planned  until  year  3  of  the  study,  based  on  current  accrual  trends  with  our 
clinical  collaborators  at  the  Boston  Medical  Center  for  different  projects,  the  number  of  target 
subjects  for  the  clinical  study  may  need  to  modified,  or  an  additional  clinical  site  may  need  to 
be  added  to  the  study  in  order  to  reach  the  n=30  target  accrual.  Additional,  we  are  currently 
considering  measuring  ex-vivo  specimens  with  hyperspectral  SFDI.  More  investigation  into 
accrual  limitations  and  other  clinical  sites  will  be  explored  over  the  next  year  to  mitigate 
changes  to  the  target  accrual. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 
There  have  been  minor  delays  in  several  of  the  subtasks,  but  we  believe  this  is  well  balanced 
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wearables”,  SPIE  Photonics  West,  San  Francisco,  California,  February  2017,  invited  oral 
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Darren  Roblyer,  Syeda  Tabassum,  Junjie  Wu,  David  J.  Waxman,  “Spatial  frequency  domain 
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-  A  new  way  to  conduct  simultaneous  multi-wavelength  (hyperspectral)  Spatial 
Frequency  Domain  Imaging  (SFDI)  has  been  developed  and  published.  This 
system  is  described  in  detail  in  section  3,  under  the  research  update  for  subtask 
3. 

■  A  new  SFDI  near  infrared  (NIR)  and  short  wave  infrared  (SWIR)  system  was 
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1.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
■  What  individuals  have  worked  on  the  project? 


Name: 

Darren  Roblyer 

Project  Role: 

PI 

Researcher  Identifier  (era  Commons 

ID): 

droblyer 

Nearest  person  month  worked: 

2 

Contribution  to  Project: 

Overall  project  management,  data  analysis,  mentorship 

Funding  Support: 

DOD 

Name: 

Irving  Bigio 

Project  Role: 

Collaborator 

Researcher  Identifier  (era  Commons 

ID): 

ijbigio 

Nearest  person  month  worked: 

1 

Contribution  to  Project: 

Mentor  and  collaborator.  Provides  technical  support  and 
feedback  for  optical  instrumentation. 

Funding  Support: 

DOD 

Name: 

David  Waxman 

Project  Role: 

Collaborator 

Researcher  Identifier  (era  Commons  ID): 

David_Waxman 

Nearest  person  month  worked: 

1 

Contribution  to  Project: 

Mentor  and  collaborator  for  small  animal  studies  and  IHC. 

Funding  Support: 

DOD 

■  Has  there  been  a  change  in  the  active  other  support  of  the  PD/PI(s)  or  senior/key  personnel  since 
the  last  reporting  period? 

■  Nothing  to  Report. 

■  What  other  organizations  were  involved  as  partners? 

■  Organization  Name:  Modulated  Imaging  Inc. 

■  Location  of  Organization:  Irvine,  CA 

■  Partner’s  contribution  to  the  project 

■  Collaboration  Amman  Mazhar  and  David  Cuccia  at  Modulated  Imaging  Inc.  provide  technical 
support  and  feedback  on  imaging  techniques  throughout  the  project. 

8.  SPECIAL  REPORTING  REQUIREMENTS 

■  Not  applicable 

9.  APPENDICES: 


p  Q  a  OOOOOOOC-OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOODOOOOOOOOOOOOOOOOO 

PoOPaooPoPooooooPoooaopoPooooooaoe>oOdooopooo*ooooPoooooooooPooOo<iooooooooooPoooGioooaooooPo 


BiomedicalOptics.SPIEDigitalLibrary.org 


Ultrafast  wavelength  multiplexed 
broad  bandwidth  digital  diffuse 
optical  spectroscopy  for  in  vivo 
extraction  of  tissue  optical  properties 


Alyssa  Torjesen 
Raeef  Istfan 
Darren  Roblyer 


SPIE. 


Alyssa  Torjesen,  Raeef  Istfan,  Darren  Roblyer,  “Ultrafast  wavelength  multiplexed  broad  bandwidth  digital 
diffuse  optical  spectroscopy  for  in  vivo  extraction  of  tissue  optical  properties,”  J.  Biomed.  Opt. 

22(3),  036009  (2017),  doi:  10.1117/1.JB0.22.3.036009. 


Downloaded  From:  http://biomedicaloptics.spiedigitallibrary.org/  on  03/10/2017  Terms  of  Use:  http://spiedigitallibrary.org/ss/termsofuse.aspx 


Journal  of  Biomedical  Optics  22(3),  036009  (March  2017) 


Ultrafast  wavelength  multiplexed  broad  bandwidth 
digital  diffuse  optical  spectroscopy  for  in  vivo 
extraction  of  tissue  optical  properties 


Alyssa  Torjesen,  Raeef  Istfan,  and  Darren  Roblyer* 
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Abstract.  Frequency-domain  diffuse  optical  spectroscopy  (FD-DOS)  utilizes  intensity-modulated  light  to  char¬ 
acterize  optical  scattering  and  absorption  in  thick  tissue.  Previous  FD-DOS  systems  have  been  limited  by  large 
device  footprints,  complex  electronics,  high  costs,  and  limited  acquisition  speeds,  all  of  which  complicate  access 
to  patients  in  the  clinical  setting.  We  have  developed  a  new  digital  DOS  (dDOS)  system,  which  is  relatively 
compact  and  inexpensive,  allowing  for  simplified  clinical  use,  while  providing  unprecedented  measurement 
speeds.  The  dDOS  system  utilizes  hardware-integrated  custom  board-level  direct  digital  synthesizers  and 
an  analog-to-digital  converter  to  generate  frequency  sweeps  and  directly  measure  signals  utilizing  undersam¬ 
pling  at  six  wavelengths  modulated  at  discrete  frequencies  from  50  to  400  MHz.  Wavelength  multiplexing  is 
utilized  to  achieve  broadband  frequency  sweep  measurements  acquired  at  over  97  Hz.  When  compared  to 
a  gold-standard  DOS  system,  the  accuracy  of  optical  properties  recovered  with  the  dDOS  system  was  within 
5.3%  and  5.5%  for  absorption  and  reduced  scattering  coefficient  extractions,  respectively.  When  tested  in  vivo , 
the  dDOS  system  was  able  to  detect  physiological  changes  throughout  the  cardiac  cycle.  The  new  FD-dDOS 
system  is  fast,  inexpensive,  and  compact  without  compromising  measurement  quality.  ©  The  Authors.  Published  by  spie 

under  a  Creative  Commons  Attribution  3.0  Unported  License.  Distribution  or  reproduction  of  this  work  in  whole  or  in  part  requires  full  attribution  of  the 
original  publication,  including  its  DOI.  [DOI:  10.1 1 17/1  .JB0.22. 3. 036009] 
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1  Introduction 

Diffuse  optical  spectroscopy  (DOS)  has  emerged  over  the  past 
three  decades  as  an  important  noninvasive  method  to  character¬ 
ize  optical  scattering  and  absorption  and,  in  turn,  calculate 
chromophore  concentrations  in  thick  tissue.1,2  Diffuse  optical 
methodologies  can  be  broadly  categorized  into  three  distinct 
measurement  types:  time  domain,  frequency  domain  (FD), 
and  continuous  wave  (CW).  This  work  focuses  on  FD-DOS, 
which  allows  for  the  separation  of  optical  scattering  and  absorp¬ 
tion  through  the  measurement  of  amplitude  attenuation  and 
phase  delay  of  photon  density  waves  generated  by  temporally 
intensity  modulated  light  sources  at  the  tissue  surface.  Prior 
FD-DOS  systems  have  utilized  either  a  single  modulation 
frequency3-9  or  a  sweep  of  multiple  frequencies5-7  in  the  radio 
frequency  (RF)  range.  Most  single  modulation  frequency  sys¬ 
tems  utilize  either  homodyne4,7,8,10  or  heterodyne3  techniques 
to  downconvert  the  high  frequency  signal  of  interest  by  mixing 
it  with  a  reference  signal.  Broad  bandwidth  DOS  instruments, 
which  utilize  light  sources  that  are  swept  through  a  range  of  RF 
frequencies,  have  typically  utilized  a  benchtop  network  analyzer 
to  generate  and  measure  signals  up  to  the  gigahertz  range.1,11 
FD-DOS  devices  based  around  network  analyzers  have  rela¬ 
tively  large  instrument  footprints,  and  the  cost  of  fabrication 
can  range  from  approximately  $30,000  to  $60,000  (based  on 
our  internal  cost  estimates),  limiting  access  to  patients  in  the 
clinic.  Although  some  analog  systems  have  been  successfully 
miniaturized,12,13  the  difficulties  involved  in  RF  system  design 


*Address  all  correspondence  to:  Darren  Roblyer,  E-mail:  roblyer@bu.edu 


potentially  limit  the  adoption  of  these  techniques.  Furthermore, 
FD-DOS  instruments  often  have  slow  measurement  acquisition 
speeds  compared  to  their  CW  counterparts,  and  although  recent 
improvements  have  been  made  in  data  acquisition  rates  for  sin¬ 
gle  modulation  frequency  measurements,14  data  acquisition 
speeds  are  limited  for  broad  modulation  frequency  bandwidth 
FD-DOS  systems.1,10,15 

We  recently  demonstrated  the  feasibility  of  FD-DOS  mea¬ 
surements  using  all  digital  signal  generation  and  detection 
using  direct  digital  synthesizer  (DDS)  integrated  circuits  and  a 
3.6-gigasample/s  two-channel  analog-to-digital  converter 
(ADC).16  These  digital  components  reduced  system  costs,  but 
measurement  acquisition  speed  was  limited  due  to  data  transfer 
bottlenecks  between  system  components.  Furthermore,  the  high 
sampling  rate  of  the  ADC  resulted  in  large  datasets  and  lengthy 
data  transfer  times.  A  second  implementation  of  a  digital  FD- 
DOS  system  tested  the  feasibility  of  utilizing  a  slower  ADC, 
sampling  at  only  25  megasamples/s,  and  utilizing  undersam¬ 
pling  techniques.17  This  substantially  reduced  system  complex¬ 
ity  but  measurement  speed  was  still  a  major  limitation  due  to 
limited  data  transfer  rates.  Here,  we  present  a  new  FD-DOS  sys¬ 
tem  that  overcomes  these  past  limitations  by  fully  integrating 
digital  signal  synthesis  and  detection  to  rapidly  acquire  multi¬ 
wavelength,  multiplexed,  broadband  frequency  sweep  measure¬ 
ments  at  a  repetition  rate  of  up  to  97  Hz.  The  digital  DOS 
(dDOS)  system  was  validated  relative  to  a  gold- standard  net¬ 
work  analyzer  based  DOS  system  and  tested  in  vivo.  This 
new  high-speed,  broad-bandwidth,  digital,  miniaturized  DOS 
system  minimizes  analog  circuitry  while  reducing  cost  and 
device  footprint,  potentially  enabling  access  to  more  patients 
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in  the  clinic  and  expanding  the  capabilities  of  this  technology, 
especially  for  monitoring  fast  physiological  changes. 


2  System  Design  and  Processing 

2.1  Hardware  Integration 

Figure  1  shows  the  core  components  and  signal/data  flow  for  the 
new  dDOS  system.  The  hardware  components  of  the  dDOS  sys¬ 
tem  are  integrated  into  a  single  motherboard  that  both  generates 
and  measures  RF  signals  while  communicating  and  transferring 
data  to  the  host  computer.  A  development  board  (ZedBoard, 
MircoZed  Zynq-7010  SoC)  was  used  as  the  central  processing 
and  data  storage  core  for  the  system.  It  contains  an  ARM 
Cortext-A9  processor  and  an  Artix  7  FPGA,  1  GB  of  RAM  to 
store  data  read  from  the  ADC,  and  a  Linux  server  to  transfer  data 
between  the  board  and  host  computer.  The  user  sets  the  desired 
frequency  sweep  settings  in  software  and  the  Linux  server  con¬ 
figures  the  DDS  boards  accordingly.  The  two-channel  ADC 
samples  both  channels  simultaneously  and  transfers  data  to 
the  RAM  on  the  MicroZed  module  via  a  parallel  low  voltage 
differential  signal  bus.  Data  are  then  transferred  from  the 
Linux  server  to  the  host  computer  via  an  Ethernet  link  at  a 
rate  of  up  to  170  Mbps. 

2.2  Light  Sources  and  Modulation 

The  dDOS  system  utilizes  NIR  laser  diodes  with  wavelengths  of 
658,  690,  785,  808,  830,  and  850  nm  that  are  housed  in  thermo¬ 
electric  cooling  modules  (Thorlabs,  LDM9T).  All  laser  diodes 
are  coupled  to  400-ywm  core  diameter  fibers  bundled  into  a  single 
ferrule.  Each  laser  diode  is  driven  with  a  DC  current,  which 
is  mixed  with  an  RF  current  using  a  bias-tee  to  intensity- 
modulated  light.  DC  current  is  supplied  by  an  eight-channel, 
high-stability  laser  diode  controller  (ILX  Lightwave,  LDC- 
390,  Newport  Corporation).  The  dDOS  system  replaces  the 
RF  current  output  functionality  of  the  network  analyzer  by 
implementing  six  DDS  chips  (AD9910,  Analog  Devices), 
each  on  individual  daughterboards.  Each  DDS  can  output  up 


to  20  mA  of  current  modulated  up  to  400  MHz.  The  RF  output 
of  each  DDS  board  is  low-pass  filtered  at  400  MHz  and  routed  to 
a  7 : 1  directional  coupler.  The  directional  coupler  sends  the 
majority  of  the  output  power  to  the  bias-tee,  and  the  lesser  out¬ 
put  from  each  directional  coupler  is  routed  to  a  power  combiner 
(Mini-Circuits,  ZBSC-615+)  and  measured  as  a  reference  signal 
with  the  second  channel  of  the  ADC. 

The  dDOS  system  is  capable  of  utilizing  the  six  individual 
DDS  chips  to  simultaneously  modulate  a  bank  of  laser  diodes  at 
offset  frequencies.  The  data  are  then  demultiplexed  in  the  FD  in 
postprocessing.  A  typical  dDOS  measurement,  for  example, 
may  sweep  through  frequencies  between  50  and  400  MHz  in 
steps  of  1  MHz  with  the  modulation  frequency  for  each  wave¬ 
length  offset  by  10  MHz.  This  multiplexing  reduces  measure¬ 
ment  duration  to  less  than  100  ms/frequency  sweep  (in  the 
slowest  case),  allowing  for  detection  of  rapid  physiological 
changes.  For  this  work,  modulation  frequencies  between  50 
and  300  MHz  or  between  50  and  400  MHz  were  used  depending 
on  the  optical  attenuation  of  the  sample  and  data  acquisition  rate 
requirements.  The  frequency  step  size  can  be  set  between  1  and 
7  MHz.  For  wavelengths  with  modulation  sweeps  starting  above 
50  MHz,  the  modulation  frequency  was  incremented  up  to  the 
maximum  frequency  then  wrapped  back  to  50  MHz  and  incre¬ 
mented  to  the  starting  frequency.  Data  were  unwrapped  such 
that  modulation  frequencies  appear  in  increasing  order. 

2.3  Detection 

Reflectance  mode  geometry  was  used  for  all  dDOS  measure¬ 
ments,  with  the  source  fiber  bundle  and  active  area  of  an  ava¬ 
lanche  photodiode  (APD)  placed  directly  on  the  surface  of  the 
phantom  at  a  source  detector  separation  of  10  to  30  mm  depend¬ 
ing  on  the  specific  measurement.  An  APD  module  (Hamamatsu 
C5658)  with  a  0.5-mm  active  area,  2.50  x  105  V/W  photoelec¬ 
tric  sensitivity,  and  a  high  band  cutoff  of  1  GHz  was  used  for 
most  phantom  measurements.  For  in  vivo  experiments,  data 
were  acquired  with  a  3 -mm  active  area  APD. 

A  14-bit,  two-channel  ADC  (Texas  Instruments,  ADS62P49) 
with  a  2  V  peak-to-peak  full-scale  input  voltage  was  used  to 


Fig- 1  Schematic  overview  of  dDOS  system.  GUI,  graphical  user  interface;  ADC,  analog-to-digital  con¬ 
verter;  FPGA,  field  programmable  gate  array;  RAM,  rapid  access  memory;  DDS,  direct  digital  synthe¬ 
sizer;  APD,  avalanche  photodiode. 
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sample  the  signal  at  each  channel  at  250  MSPS.  The  electrical 
output  of  the  APD  was  high-pass  filtered  at  41  MHz  and  routed 
to  the  first  channel  of  the  ADC,  while  the  combined  reference 
signal  from  the  power  combiner  was  routed  to  the  second  chan¬ 
nel  of  the  ADC.  The  user  can  select  the  number  of  samples  col¬ 
lected  by  each  channel  of  the  ADC  at  each  frequency  step  (e.g., 
4096,  8192, 16384. . .  2n).  Increasing  the  number  of  samples  per 
step  substantially  increases  data  transfer  and  processing  time. 
Both  channels  of  the  ADC  run  off  of  the  same  clock  and  are 
sampled  simultaneously. 


2.4  Signal  Processing 

All  signal  processing  is  performed  in  MATLAB  (R2014b, 
Mathworks  Natick,  Massachusetts).  Since  the  two-channel 
ADC  has  a  sampling  rate  of  250  MSPS,  all  signals  modulated 
above  125  MHz  are  aliased.  Given  that  the  modulation  fre¬ 
quency  is  known  for  each  wavelength  at  each  frequency  step, 
the  signal  of  interest  is  easily  located  in  the  baseband.  The 
processing  algorithm  eliminates  modulation  frequencies  that 
are  multiples  of  the  Nyquist  frequency  and  removes  any  mea¬ 
surements  in  which  multiple  wavelengths  are  aliased  to  the  same 
frequency  in  the  baseband  (e.g.,  signals  at  120  and  130  MHz  are 
both  read  as  120  MHz  and  eliminated). 

A  fast  Fourier  transform  (FFT)  with  a  rectangular  window 
was  used  at  each  frequency  step  to  calculate  the  amplitude 
and  phase  for  the  reference  and  sample  channels.  The  appropri¬ 
ate  frequency  bin  was  located  for  each  wavelength  at  each 
modulation  frequency  step.  The  magnitude  of  the  FFT  of  the 
measurement  channel  relative  to  that  of  the  reference  channel 
was  considered  the  raw  amplitude  value,  and,  similarly,  the 
phase  offset  between  the  measurement  and  reference  channels 
was  considered  the  raw  phase  value. 


2.5  Measurement  Calibration  and  Optical  Property 
Extraction 

The  raw  amplitude  and  phase  represent  the  amplitude  attenu¬ 
ation  and  phase  delay  induced  both  by  the  tissue  and  the 
instrument.  In  order  to  remove  the  instrument  response,  a 
measurement  was  taken  on  a  silicone  phantom  with  known 
optical  properties,  where  calibration  optical  properties  were 
obtained  using  multidistance  measurements.  An  analytical  sol¬ 
ution  to  the  PI  approximation  of  the  Boltzman  transport  equa¬ 
tion  with  boundary  conditions  for  semi-infinite  geometry  was 
employed  as  a  forward  model  for  the  phantom  measurement.5 
Calibration  factors  for  amplitude  and  phase  were  determined 
by  ratiometrically  or  differentially  comparing  the  measured 
amplitude  and  phase,  respectively,  to  the  theoretical  amplitude 
and  phase  from  the  forward  model.  These  calibration  factors 
were  then  applied  to  the  subsequent  raw  amplitude  and 
phase  measurements  in  order  to  obtain  calibrated  amplitude 
and  phase. 

An  iterative  least-squares  fitting  algorithm  was  employed  to 
minimize  the  error  between  the  forward  model  calculated  and 
measured  calibrated  amplitude  and  phase  to  extract  absorption 
and  scattering  (//s')  parameters  at  each  wavelength.5 
Chromophore  extinction  coefficients  were  obtained  from 
Zijlstra  and  Buursma.18  The  Beer-Lambert  law  was  used  to 
extract  chromophore  concentrations  (oxy-  and  deoxyhemoglo¬ 
bin)  from  multiwavelength  yt/a  data.13 


3  System  Characterization 

3.1  System  Overview 

The  completed  dDOS  system  is  shown  in  Fig.  2,  where  the 
motherboard,  DDS  daughterboards,  MicroZed  module,  ADC, 
and  RF  output  connections  are  pictured.  The  primary  hardware 
is  housed  in  an  enclosure  measuring  10  X  12x6  in. 

3.2  Signal  to  Noise 

Signal  to  noise  was  measured  on  a  breast-simulating  silicone 
optical  phantom.  SNR  was  compared  for  measurements  taken 
using  either  a  single  laser  or  all  six  lasers  modulated  simultane¬ 
ously.  A  comparison  was  also  done  for  measurements  taken 
using  either  4096  or  8192  samples.  The  noise  floor  was  calcu¬ 
lated  as  the  average  noise  level,  in  dBc,  relative  to  the  frequency 
of  interest.  Figure  3  shows  that  when  all  six  lasers  were  modu¬ 
lated  simultaneously,  the  noise  floor  was  -46.3  dBc,  whereas 
utilizing  a  single  laser  yielded  a  lower  noise  floor  of  —50. 1  dBc. 
When  increasing  the  sample  length  to  8192  samples  per  fre¬ 
quency  step,  the  noise  floor  was  reduced  to  -49.2  and 
-54.5  dBc  for  multiplexed  and  single  wavelength  measure¬ 
ments,  respectively.  Although  the  noise  floor  was  lower  with 
8192  samples  per  frequency  step,  acquisition  and  data  transfer 
times  of  these  measurements  were  substantially  slower  than  of 
those  with  4096  samples  per  frequency  step,  and  there  was  no 
improvement  in  accuracy  of  optical  properties  (e.g.,  the  accu¬ 
racy  of  10  measurements  on  two  phantoms  with  4096  samples 
per  frequency  step  was  1.29%  compared  to  1.85%  with  8192 
samples  per  frequency  step).  4096  samples  per  frequency 
step  were  collected  for  all  subsequent  measurements  to  optimize 
speed. 

3.3  Accuracy  of  Optical  Property  Extractions 

An  example  of  calibrated  (black)  and  fit  (red)  amplitude  and 
phase  measured  from  50  to  400  MHz  in  steps  of  1  MHz  at 
each  wavelength  is  shown  in  Fig.  4.  Data  were  calibrated 
using  a  silicone  phantom  with  known  optical  properties  in 
order  to  remove  the  instrument  response.  Measurements  at 
658  and  850  nm  have  somewhat  higher  noise  levels  compared 
to  other  wavelengths,  likely  due  to  the  lower  modulation  index 
achieved  for  these  laser  diodes. 


Fig.  2  dDOS  core  electronics. 
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Fig.  3  Noise  floor  when  collecting  4096  samples  per  frequency  step  in  (b)  a  single  wavelength  meas¬ 
urement  or  (a)  a  multiplexed  measurement  with  a  10  MHz  offset  between  wavelengths. 
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each  system  with  a  15 -mm  source-detector  separation  and  all 
six  wavelengths  modulated  at  frequencies  between  50  and 
400  MHz  with  1-MHz  frequency  steps.  Measurements  on 
one  highly  attenuating  phantom  at  658,  690,  and  850  nm 
were  excluded  due  to  phase  errors  resulting  from  low  detected 
signal  levels.  The  mean  of  the  absolute  difference  in  optical 
properties  for  all  measured  phantoms  between  the  two  systems 
was  5.3%  and  5.5%  for  and  yt/',  respectively;  when  frequency 
sweeps  from  50  to  300  MHz  were  used,  the  mean  accuracy  in 
optical  properties  measured  with  the  two  systems  was  4.9%  and 
6.0%  for  and  //s',  respectively. 

Figure  5  shows  scatterplots  (top  panels)  and  Bland-Altman 
plots  (bottom  panels)  for  fid  and  //s'  measured  with  the  dDOS 
and  benchtop  DOS  systems.  For  the  Bland-Altman  plots,  the 
differences  between  optical  properties  with  the  dDOS  and  bench- 
top  system  are  plotted  on  the  y-axis  versus  average  optical  proper¬ 
ties  for  each  system  on  the  x-axis,  for  each  wavelength  and 
phantom.  96%  of  fid  measurements  and  98%  of  {i[  measurements 
fell  within  1 .96  standard  deviations  of  the  mean  difference,  indi¬ 
cating  acceptable  equivalence  between  the  two  systems.19  \id  and 
//s'  both  tend  to  skew  slightly  lower  for  the  dDOS  system  com¬ 
pared  to  the  benchtop  DOS  system.  pid  values  had  the  best  agree¬ 
ment  between  systems  when  the  mean  absorption  coefficient  was 
<0.015  mm-1,  while  //s'  had  the  best  agreement  between  systems 
when  the  mean  reduced  scattering  coefficient  was  >1  mm-1.  In 
healthy  breast  tissue,  fia  in  the  NIR  is  typically  <0.01  mm-1  and 
jil  in  the  NIR  is  typically  >0.8  mm-1.20 
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Fig.  4  Amplitude  (in  arbitrary  units,  A.U.)  and  phase  measured  with 
dDOS  system  at  each  wavelength.  Red  lines  are  calibrated  data  and 
black  lines  are  model  fits. 


dDOS  optical  property  extraction  accuracy  was  determined 
by  measuring  a  set  of  10  tissue-simulating  silicone  phantoms 
with  varying  amounts  of  titanium  dioxide  and  nigrosin,  as 
the  scattering  and  absorbing  agents,  respectively,  using  both 
the  dDOS  and  the  network  analyzer-based  benchtop  DOS  sys¬ 
tems.  For  each  phantom,  10  FD  measurements  were  taken  with 


3.4  Precision 

The  stability  of  amplitude,  phase,  and  optical  properties  was 
characterized  during  a  2-h  drift  test.  Measurements  were 
taken  at  ~l-min  intervals  on  a  silicone  optical  phantom  using 
658,  690,  785,  808,  830,  and  850  nm  lasers  modulated  simulta¬ 
neously  in  frequency  sweeps  from  50  to  400  MHz  with  a  15 -mm 
source-detector  separation.  Amplitude  and  phase  were  exam¬ 
ined  for  the  duration  of  the  drift  test  for  the  850-nm  laser  modu¬ 
lated  at  50  MHz  and  were  plotted  in  Fig.  6.  For  this  wavelength 
and  modulation  frequency,  the  standard  deviation  of  the  phase 
was  0.38  deg  and  the  standard  deviation  of  the  amplitude  rela¬ 
tive  to  the  mean  amplitude  was  0.81%. 

System  stability  was  further  analyzed  for  all  wavelengths  and 
modulation  frequencies  used  in  the  drift  test,  and  results  are 
shown  in  Table  1.  Amplitude  and  optical  property  stability 
were  defined  as  the  standard  deviation  divided  by  the  mean 
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Fig.  5  (a  and  b)  Scatterplots  and  (c  and  d)  Bland-Altman  plots  of  optical  properties  obtained  with  the 
dDOS  system  compared  to  those  obtained  with  the  benchtop  system.  Top  and  bottom  dashed  lines 
represent  1.96  standard  deviations  from  the  mean  difference;  center  dashed  lines  represent  the 
mean  difference. 


Amplitude  Drift  at  50  MHz  and  850  nm 
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Fig.  6  2-h  drift  test  showing  dDOS  system  stability.  Measurements 
were  taken  approximately  every  minute  and  show  variation  from 
value  at  the  first  time  point  ( t  =  0).  (a)  Amplitude  drift  at  50  MHz 
and  850  nm  and  (b)  phase  drift  at  50  MHz  and  850  nm. 


Table  1  System  stability  during  a  2-h  drift  test  with  measurements 
taken  approximately  every  minute. 


Average  stability  across  all  wavelengths  and  modulation  frequencies 

Amplitude 

1 .0% 

Phase 

1.7  deg 

f^a 

2.0% 

0.8% 

at  each  modulation  frequency,  averaged  for  all  wavelengths. 
Phase  stability  was  defined  as  the  standard  deviation  throughout 
the  drift  test  averaged  for  all  wavelengths  and  modulation 
frequencies.  The  stability  in  fia  and  //s'  was  2.0%  and  0.8%, 
respectively. 


3.5  Speed 

Most  prior  FD-DOS  instruments  performed  frequency  sweeps 
with  each  wavelength  modulated  individually  in  a  sequential 
manner.5,17  The  dDOS  system  utilizes  wavelength  multiplexing, 
which  eliminates  switching  time  between  lasers  and  reduces 
data  acquisition  time,  resulting  in  the  ability  to  measure  fast 
physiological  changes.  The  rate  at  which  frequency  sweep  mea¬ 
surements  can  be  performed  is  dependent  on  the  range  of  modu¬ 
lation  frequencies,  the  step  size,  and  the  number  of  samples 
acquired  at  each  frequency  step.  When  multiplexing  six  laser 
wavelengths,  sweeping  from  50  to  300  MHz  with  7  MHz 
steps,  and  acquiring  4096  samples  per  frequency  step,  the  maxi¬ 
mum  speed  achieved  for  dDOS  system  was  97.2  Hz.  Increasing 
the  frequency  sweep  range  from  50  to  400  MHz  reduces  the 
measurement  rate  to  68.1  Hz.  When  decreasing  the  frequency 
step  size  to  1  MHz,  the  50  to  300  MHz  sweep  can  be  acquired 
at  13.9  Hz,  and  the  50  to  400  MHz  sweep  can  be  acquired  at 
9.8  Hz.  These  measurement  speeds  allow  for  FD-dDOS  data 
acquisition  that  can  capture  fast  physiological  changes,  includ¬ 
ing  during  the  cardiac  cycle. 

To  test  the  maximum  measurement  speed  attainable  with  the 
dDOS  system,  measurements  were  also  taken  with  a  single 
modulation  frequency.  When  modulating  at  only  50  MHz,  a 
measurement  rate  of  2557  Hz  can  be  achieved.  Although  this 
is  not  the  preferred  data  acquisition  mode,  because  utilizing 
additional  modulation  frequencies  improves  model  fit  accuracy, 
it  provides  a  useful  comparison  to  other  FD  systems,  which  uti¬ 
lize  single  modulation  frequencies.3,4,8,10,14,21 

A  summary  of  the  performance  specifications  is  given  in 
Table  2. 
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Table  2  dDOS  performance  specifications. 


dDOSI  performance  summary 

SNR  (typical) 

46.3  dBc 

fia  accuracy 

5.3% 

Us  accuracy 

5.5% 

Phase  stability 

1.7  deg 

Amplitude  stability 

1 .0% 

Measurement  rate  (50  to  300  MHz  sweeps) 

97.2  Hz 

Maximum  measurement  rate  (single  frequency) 

2557  Hz 

4  In  Vivo  Testing 

All  in  vivo  human  measurements  were  carried  out  under  a  pro¬ 
tocol  approved  by  the  Boston  University  Institutional 
Review  Board. 

4.1  Cuff  Occlusion  Measurement 

An  in  vivo  cuff  occlusion  measurement  was  performed  using  a 
blood  pressure  cuff  placed  on  the  upper  arm  of  a  male  volunteer, 
while  optical  measurements  were  taken  on  the  forearm.  The  cuff 
was  loosely  secured  in  order  to  measure  baseline  molar  concen¬ 
trations  of  oxyhemoglobin  (ct02Hb)  and  deoxyhemoglobin 
(ctHHb).  The  cuff  was  then  was  inflated  to  200  mmHg  for 
~2  min,  and  then  released  and  data  were  acquired  for  ~5  addi¬ 
tional  minutes.  Frequency  sweeps  between  50  and  400  MHz, 
with  frequency  steps  of  1  MHz,  were  used  to  modulate  all 
six  laser  diodes  simultaneously.  Measurements  were  taken 
approximately  every  8  s  for  10  min  and  a  source-detector  sep¬ 
aration  of  10  mm  was  used  for  each  measurement. 

Figure  7  shows  ct02Hb  and  ctHHb  throughout  the  cuff 
occlusion  test.  At  baseline  average,  ct02Hb  and  ctHHb  were 

41.6  and  41.2  piM,  respectively.  At  the  time  of  cuff  release, 
ct02Hb  had  dropped  to  19.3  //M  and  ctHHb  had  increased  to 

87.6  //M.  After  cuff  release,  ct02Hb  rebounded  to  a  local  maxi¬ 
mum  of  67.7  //M  and  ctHHb  dropped  to  36.5  fiM,  before  both 
chromophore  concentrations  began  to  return  toward  their 


Cuff  Occlusion  Test 


Fig.  7  Forearm  cuff  occlusion  test.  Cuff  was  occluded  during  time 
interval  between  dashed  lines.  Measurements  were  taken  approxi¬ 
mately  every  8  s  and  a  three-point  moving  average  is  plotted  for 
each  chromophore. 


baseline  values.  These  results  are  consistent  with  reduced 
arterial  supply  and  vascular  drainage  during  cuff  occlusion,  fol¬ 
lowed  by  a  rush  of  arterial  blood  into  the  forearm  after  cuff 
release. 

4.2  Rapid  In  Vivo  Measurements 

To  test  the  ability  of  the  dDOS  system  to  measure  physiological 
changes  at  the  cardiac  rate,  frequency  sweep  measurements 
were  taken  on  the  finger  of  a  male  volunteer.  Measurements 
were  taken  in  reflectance  mode  on  the  index  finger  with  a 
10-mm  source-detector  separation.  Frequency  sweeps  from 
50  to  300  MHz  with  7  MHz  steps  were  used  to  modulate  all 
six  laser  diodes  simultaneously.  About  350  frequency  sweeps 
were  taken  with  a  measurement  rate  of  97.2  Hz. 

An  electrocardiogram  (ECG)  signal  was  measured  simulta¬ 
neously  with  the  optical  signal  and  synced  using  an  external  data 
acquisition  board.  Figure  8  shows  five-point  moving  averages  of 
//a,  //s',  ctTHb,  ct02Hb,  and  ctHHb  from  the  high-speed  finger 
measurement  temporally  aligned  with  the  ECG.  //a  values  fluc¬ 
tuate  temporally  with  the  ECG  signal,  likely  due  to  increased 
blood  volume  in  the  finger  after  ventricular  contraction. 
Similarly,  ctTHb  and  ct02Hb  vary  temporally  with  the  ECG  sig¬ 
nal,  while  ctHHb  shows  minimal  fluctuation  throughout  the 
cardiac  cycle.  Vertical  dashed  lines  mark  the  peak  of  the 
ECG  R-wave,  the  start  of  the  upstroke  of  the  pulsatile  waveform, 
and  the  dicrotic  notch.  There  is  a  155-ms  lag  in  the  upstroke  of 
the  //a  signal  after  the  R-wave  of  the  ECG,  due  to  the  transmis¬ 
sion  time  of  the  pressure  wave  to  the  peripheral  vasculature.  The 
delay  is  often  referred  to  as  the  pulse  transit  time  and  previously 
shown  to  range  between  200  and  300  ms  for  photoplethysmo- 
graphic  signals  measured  in  the  fingertip.22  The  dicrotic  notch, 
which  is  correlated  with  end  of  the  systolic  ejection  period,  is 
likely  correlated  with  the  dip  in  //a,  ctTHb  and  ct02Hb  marked 
by  the  rightmost  vertical  line.  In  this  case,  the  dichrotic  notch 
appears  roughly  303  ms  after  the  upstroke  of  the  pulsatile  wave¬ 
form,  in  agreement  with  other  report  values  in  the  literature.23,24 
Fluctuations  in  //a,  ctTHb,  and  ct02Hb  appear  to  occur  at  the 
same  frequency  as  the  ECG  signal,  suggesting  that  the  dDOS 
system  is  capable  of  detecting  pulsatile,  oxygenated,  arterial 
blood  flow. 

4.3  Chromophore  Maps  in  Healthy  Breast  Tissue 

Optical  properties  of  breast  tissue  were  measured  on  a  healthy 
female  volunteer.  Measurements  were  taken  with  the  dDOS  sys¬ 
tem  using  658,  690,  785,  808,  830,  and  850  nm  lasers  with  a  28- 
mm  source-detector  separation  and  modulation  frequencies 
from  50  to  400  MHz.  36  point  measurements  were  taken  in 
which  the  source-detector  pair  was  manually  scanned  across 
the  right  breast  in  a  6  X  6  cm  grid  pattern  with  1  cm  spacing 
between  measurement  points.  A  quadrant  of  the  areolar  complex 
was  located  at  the  top  left  comer  of  the  imaging  grid.  Optical 
properties  were  calculated  and  chromophore  values  determined 
for  each  data  point.  Data  from  the  658-  and  850-nm  wavelengths 
were  excluded  due  to  noise  in  the  phase  measurement  at  these 
wavelengths.  Figure  9  shows  a  chromophore  map  for  ctTHb 
with  linear  interpolation  between  measurement  points.  The 
chromophore  map  was  superimposed  on  a  3-D  scanned  bust 
model  in  order  to  show  the  approximate  relative  position  of 
the  measurement.  There  is  a  region  of  increased  ctTHb  at  the 
top  left  comer  of  the  chromophore  map,  likely  due  to  increased 
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Fig.  8  High-speed  dDOS  measurements  were  taken  on  the  index  fin¬ 
ger.  An  ECG  signal  was  obtained  simultaneously  with  the  optical 
measurements  and  temporally  aligned  with  optical  property  and 
chromophore  values.  From  top  to  bottom,  plots  represent  the  ECG, 
lia,  /4,  ctTHb,  ct02Hb,  and  ctHHb.  The  leftmost  vertical  dashed 
line  represents  the  approximate  time  of  the  peak  of  the  ECG  IT- 
wave,  the  middle  dashed  line  represents  the  approximate  time  of 
the  start  of  the  upstroke  of  the  optical  signal  following  the  arterial 
transit  time,  and  the  rightmost  dashed  line  represents  the  approximate 
timing  of  the  dicrotic  notch  based  on  visual  inspection. 


blood  supply  to  the  areola  as  compared  to  standard  breast  tissue, 
which  has  been  observed  in  prior  DOS  studies.25 

5  Discussion 

A  new  dDOS  system  was  developed  to  increase  access  to 
patients,  while  adding  new  capabilities  in  measurement  speed. 
Traditional  network  analyzer-based  FD-DOS  systems  are  large 
and  costly,  and  miniaturization  of  analog  RF  circuits  is  techni¬ 
cally  challenging.  The  hardware-integrated  dDOS  system  devel¬ 
oped  here  utilizes  direct  digital  signal  generation  and  detection 
and  wavelength  multiplexing,  with  frequency  sweep  measure¬ 
ments  rates  up  to  97  Hz,  allowing  rapid  physiological  changes 
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Fig.  9  ctTHb  measured  on  the  right  breast  of  a  healthy  volunteer.  A 
quadrant  of  the  nipple  is  located  in  the  lower  left  corner  of  the  image. 


to  be  detected  well  above  the  cardiac  rate.  Previous  systems 
achieving  high  sampling  rates  have  utilized  a  small  number 
of  discrete  modulation  frequencies  and  wavelengths,14,21 
whereas  the  dDOS  system  sweeps  a  range  of  hundreds  of  modu¬ 
lation  frequencies  and  utilizes  six  discrete  wavelengths. 

The  dDOS  system  performance  is  comparable  to  previous 
systems.  For  example,  Pham  et  al.  reported  an  accuracy  of  ±5% 
and  ±3%  for  / ia  and  yw',  respectively,  whereas  the  dDOS  system 
had  an  accuracy  of  ±5.3%  and  ±5.5%  for  the  same  parameters. 
The  dDOS  design  has  similarities  to  that  of  the  modular  digital 
FD  diffuse  optical  monitor  developed  by  Durduran  et  al.,  which 
utilized  digital  signal  synthesis  and  a  dual-channel  ACD  to  sam¬ 
ple  and  downconvert  the  signal  from  three  laser  wavelengths  at  a 
single  modulation  frequency;21  however,  our  design  differs  in 
that  we  use  six  multiplexed  laser  wavelengths  modulated  at  hun¬ 
dreds  of  discrete  frequencies  in  the  RF  range.  Zimmerman  et  al. 
recently  developed  a  FD  NIR  high  speed  imaging  system  that 
utilized  two  wavelengths,  two  modulation  frequencies,  and  an 
array  of  detectors;  with  this  system,  a  measurement  rate  of 
90  Hz  was  achieved,  compared  to  the  97-Hz  measurement 
rate  achieved  in  the  dDOS  system  with  six  wavelengths  and 
a  sweep  of  modulation  frequencies. 14 

In  the  future,  the  portability  of  the  system  can  be  further 
improved  by  implementing  miniaturized  fiber-coupled  laser 
diodes  and  miniaturized  current  driver  modules  that  can  be 
housed  in  a  portable  device  enclosure.  This  will  help  to  more 
easily  measure  patients  in  a  clinical  setting  for  a  range  of  appli¬ 
cations,  such  as  chemotherapy  monitoring  or  rapid  hemo¬ 
dynamic  analysis.  The  speed  of  the  system  also  opens  the 
possibility  of  monitoring  rapid  physiological  changes  at  the  car¬ 
diac  rate,  which  may  be  useful  in  monitoring  hemodynamics 
related  to  conditions,  such  as  peripheral  artery  disease. 


6  Conclusions 

We  have  developed  an  inexpensive  clinical  FD-DOS  platform 
based  on  DDSs  and  ADCs  implemented  on  custom  hard- 
ware-integrated  PCBs.  The  system  has  high  stability  and  accu¬ 
racy  in  phantom  measurements  and  is  capable  of  detecting  rapid 
physiological  fluctuations  above  the  cardiac  rate  using  six  wave¬ 
lengths  and  a  sweep  of  modulation  frequencies.  The  dDOS 
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system  may  improve  access  to  patients  while  helping  to  open 
new  areas  of  investigation. 
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Abstract.  Spatial  frequency  domain  imaging  (SFDI)  is  a  wide-field  diffuse  optical  technique  used  to  obtain  opti¬ 
cal  properties  and  chromophore  concentrations  in  highly  scattering  media,  such  as  biological  tissue.  Here,  we 
present  a  method  for  rapidly  acquiring  multispectral  SFDI  data  by  modulating  each  illumination  wavelength  at  a 
different  temporal  frequency.  In  the  remitted  signal,  each  wavelength  is  temporally  demodulated  and  processed 
using  conventional  SFDI  techniques.  We  demonstrate  a  proof-of-concept  system  capable  of  acquiring  wide-field 
maps  (2048  x  1 536  pixels,  8.5  x  6.4  cm)  of  optical  properties  at  three  wavelengths  in  under  2.5  s.  Data  acquired 
by  this  method  show  a  good  agreement  with  a  commercial  SFDI  imaging  system  (with  an  average  error  of  1 3%  in 
absorption  and  8%  in  scattering).  Additionally,  we  show  that  this  strategy  is  insensitive  to  ambient  lighting  con¬ 
ditions,  making  it  more  practical  for  clinical  translation.  In  the  future,  this  technique  could  be  expanded  to  tens  or 
hundreds  of  wavelengths  without  increasing  acquisition  time.  ©  The  Authors.  Published  by  SPIE  under  a  Creative  Commons 
Attribution  3.0  Unported  License.  Distribution  or  reproduction  of  this  work  in  whole  or  in  part  requires  full  attribution  of  the  original  publication,  including 
its  DOI.  [DOI:  10.1 1 17/1  .JBO.22.7.076019] 
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1  Introduction 

At  near-infrared  wavelengths,  biological  tissue  is  highly  scatter¬ 
ing,  making  high-resolution  images  of  tissue  difficult  at  depths 
greater  than  about  1  mm.1  However,  multiply  scattered  photons 
carry  information  about  the  absorption  coefficient  (pi.d)  and  the 
reduced  scattering  coefficient  (//s'),  which  can  reveal  the  function 
and  structure  of  the  sample.  From  these  optical  properties,  the 
concentration  of  biologically  relevant  chromophores,  such  as 
oxy-  and  deoxyhemoglobin,  lipid,  and  water,  can  be  deduced. 
Knowledge  of  these  chromophore  concentrations  provides  a 
wealth  of  information  about  the  composition  and  metabolic 
activity  of  tissue,  which  has  proven  useful  in  oncology,2,3 
small  animal  imaging,4  exercise  physiology,5  and  imaging  the 
brain.6,7  Most  diffuse  optics  techniques  use  point-by -point  scan¬ 
ning  to  build  up  a  map  of  tissue  optical  properties,  but  this  strat¬ 
egy  limits  the  area  and  resolution  with  which  tissue  can  be 
probed  in  a  clinically  relevant  period  of  time.8  Efforts  to  increase 
speed  and  resolution  of  point-by-point  systems  typically  involve 
increasing  the  number  of  sources  and  detectors  (often  into  the 
hundreds),  resulting  in  devices  that  are  expensive,  bulky,  and 
specific  to  particular  anatomical  sites.9,10  Recently,  a  new  tech¬ 
nique  called  spatial  frequency  domain  imaging  (SFDI)  has  been 
shown  to  provide  wide-field  determinations  of  optical  properties 
of  scattering  media  without  point-by-point  scanning.11,12  SFDI 
involves  the  illumination  of  tissue  with  spatially  modulated  light 
and  the  recording  of  diffuse  reflectance  remitted  by  the  tissue 
using  a  camera.  The  response  of  the  tissue  in  the  spatial  fre¬ 
quency  domain  can  be  related  to  the  spatial  domain  by  an 
inverse  Fourier  transform,  resulting  in  2-D  maps  of  the  diffuse 
reflectance.  The  resulting  reflectance  maps  can  then  be  used  to 
determine  optical  properties  at  each  pixel  of  the  image. 
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The  first  described  SFDI  processing  technique  used  a  mini¬ 
mum  of  six  image  acquisitions  per  wavelength.  At  each  of  at 
least  two  spatial  frequencies,  images  of  the  sample  were 
acquired  with  three  spatially  offset  phases,  which  were  demodu¬ 
lated  to  obtain  the  envelope  of  the  diffusely  reflected  light.12 
This  envelope  contains  information  about  the  structure  and  com¬ 
position  of  the  tissue  under  investigation.  During  demodulation, 
DC  background  illumination  is  subtracted  from  the  resulting 
images,  rendering  this  technique  theoretically  insensitive  to 
ambient  light  so  long  as  it  remains  constant  during  the  acquis¬ 
ition.  However,  the  collection  of  multiple  images,  as  well  as 
the  time  needed  to  switch  between  wavelengths,  results  in 
acquisitions  times  that  typically  stretch  into  the  tens  of  seconds. 
Recently,  more  advanced  demodulation  techniques  have 
reduced  the  number  of  images  needed  to  two  per  wavelength,13 
or  one  per  wavelength  in  a  technique  termed  “single-shot  optical 
properties”  (SSOP). 14,15  This  method  exploits  the  fact  that  sinus¬ 
oidal  patterns  of  light  contain  at  least  two  spatial  frequencies:  the 
DC  offset  with  a  spatial  frequency  of  0  and  the  AC  component  at 
the  projected  spatial  frequency.  Low-pass  filtering  of  the  image 
effectively  separates  out  the  DC  component,  while  the  envelope 
of  the  higher  frequency  component  can  be  computed  using 
a  Hilbert  transform,  or  other  single-sideband  demodulation 
technique.16  The  single-shot  processing  method  reduces  the 
number  of  images  that  need  to  be  acquired  at  each  wavelength 
from  six  to  one,  but  still  requires  that  each  wavelength  be 
collected  sequentially.  Furthermore,  this  technique  is  more 
sensitive  to  changes  in  ambient  light  than  the  original  method, 
posing  a  challenge  for  fast-paced  clinical  environments,  where 
lighting  conditions  may  be  constantly  changing. 

Here,  we  present  a  new  SFDI  acquisition  method  that  enables 
the  rapid  collection  of  optical  properties  simultaneously  from  a 
large  number  of  wavelengths  while  being  insensitive  to  ambient 
light.  The  method  borrows  from  a  technique  previously  used  to 
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enhance  hyperspectral  fluorescence  imaging  speed  by  modulat¬ 
ing  each  excitation  wavelength  at  different  temporal  rates.17  In 
this  study,  we  use  three  wavelengths  of  light  that  are  temporally 
modulated  at  different  frequencies  and  then  used  to  project  a 
single  sinusoidal  illumination  pattern.  We  call  this  method  “tem¬ 
porally  modulated  SFDI”  (TM-SFDI).  By  collecting  a  sequence 
of  images  of  the  sample  while  it  is  illuminated  with  time-varying 
light,  each  pixel  of  the  video  can  be  temporally  demodulated  by 
determining  the  strength  of  the  reflected  signal  at  each  modu¬ 
lation  frequency.  The  result  of  temporal  demodulation  is  a  single 
spatially  modulated  image  at  each  wavelength,  which  can  then 
be  processed  using  the  SSOP  method.  Because  this  technique  is 
only  sensitive  to  light  modulated  at  a  specific  frequency,  it  auto¬ 
matically  rejects  ambient  light,  unlike  conventional  single-shot 
processing.  Additionally,  because  there  is  no  need  to  switch 
between  wavelengths  or  projection  patterns,  acquisition  speeds 
are  only  limited  by  the  frame  rate  of  the  camera  and  the  number 
of  video  frames  needed  for  successful  temporal  demodulation. 
In  Sec.  2,  we  will  describe  how  the  TM-SFDI  system  was 
realized;  Sec.  3  demonstrates  that  it  is  capable  of  accurately 
determining  optical  properties,  is  insensitive  to  ambient  light, 
and  can  be  used  to  image  complex  samples. 

2  Materials  and  Methods 

2.1  Hardware 

For  this  proof-of-concept  system,  three  wavelengths  of  light 
were  used:  519,  652,  and  740  nm.  One  high-power  LED  per 
wavelength  (Cree  XLamp,  Cree,  Inc.  Durham,  North  Carolina) 
was  driven  at  700  mA  (652  and  740  nm)  or  1000  mA  (519  nm) 
with  independent  constant  current  sources.  Light  from  the 
LEDs  was  modulated  in  time  using  a  pulse-width  modulation 
sequence  controlled  by  a  microcontroller  (Arduino  UNO). 
Modulation  frequencies  were  set  at  5.21,  8.68,  and  12.15  Hz, 
each  of  which  falls  into  a  single  bin  of  a  32  point  discrete 
Fourier  transform  when  a  sampling  period  of  18  ms 


(55.6  frames/s)  is  used.  Frequencies  occupying  a  single  bin 
were  important  for  reducing  cross-talk  between  wavelengths. 
These  frequencies  were  also  chosen  to  avoid  harmonic  content 
from  lower  frequencies  influencing  higher  frequency  modula¬ 
tion.  Light  from  the  LEDs  was  roughly  collimated  by  a  lens 
over  each  emitter  and  sent  through  a  transparency  printed 
with  a  sine  wave  pattern  using  a  laser  jet  printer  [Fig.  1(a)]. 
The  spatiotemporally  modulated  light  was  directed  onto  a  sam¬ 
ple,  and  the  reflected  light  was  recorded  by  a  16  bit  scientific 
CCD  camera  (Zyla  4.2,  Andor,  Inc.)  with  an  8.5  x  6.4  cm  field- 
of-view  (FoV).  Crossed  polarizers  were  used  to  minimize  the 
effect  of  specular  reflection.  Modulation  and  acquisition  timings 
were  controlled  by  a  custom  Arduino  program. 

2.2  Processing 

Videos  of  temporally  modulated  light  reflected  from  the  sample 
were  acquired  at  55.6  frames/s  for  2.3  s  (128  frames)  for  offline 
processing.  A  Fourier  transform  in  time  of  each  pixel  of 
the  image  was  performed  to  separate  the  wavelengths.  The 
amplitude  of  the  spectrum  at  each  modulation  frequency 
corresponded  to  the  intensity  of  the  light  reflected  by  that 
specific  wavelength  at  that  specific  spatial  location  [Fig.  1(b)]. 
Recording  the  strength  of  the  peaks  at  each  pixel  effectively  sep¬ 
arates  each  wavelength,  while  simultaneously  rejecting  back¬ 
ground  light  (Fig.  2).  Once  an  image  of  each  wavelength  was 
extracted,  the  images  were  rotated  to  ensure  that  the  stripes  were 
parallel  to  the  top  and  bottom  of  the  frame,  and  the  optical  prop¬ 
erties  were  determined  by  previously  published  techniques.14 
Briefly,  each  image  was  split  into  a  DC  and  AC  component 
by  taking  a  spatial  Fourier  transform  of  each  line  of  the 
image  and  filtering  in  the  spatial  frequency  domain.  The  DC 
image  was  generated  by  removing  a  band  of  frequencies  sur¬ 
rounding  the  projected  AC  frequency,  leaving  the  high  and 
low  frequency  content  intact.  The  inverse  Fourier  transform 
of  this  spectrum  was  used  to  determine  the  diffuse  reflectance 
of  the  sample  under  planar  illumination.  The  AC  image  was 


(a)  |  Sample 


(b) 


Frequency 


Fig.  1  (a)  System  diagram.  Three  LEDs  were  driven  by  independent  constant  current  sources  (D1 ,  D2, 
D3)  and  were  temporally  modulated  by  a  microcontroller.  The  light  was  roughly  collimated  by  three  lenses 
(CL)  and  polarized  (PI )  before  passing  through  a  patterned  transparency  (T).  Before  being  imaged  by  the 
camera,  reflected  light  passed  through  a  second  polarizer  (P2)  set  at  90  deg  to  PI  to  reject  specular 
reflection.  Camera  FoV  is  indicated  by  dashed  lines,  (b)  Schematic  of  data  processing.  A  video  was 
captured  from  the  camera  and  the  variation  of  each  pixel  in  time  was  plotted.  An  FFT  of  this  time  series 
was  taken  and  the  strength  of  the  signal  at  the  frequencies  corresponding  to  the  modulation  of  the  LEDs 
was  used  to  build  up  a  2-D  image  for  each  wavelength.  These  temporally  demodulated  images  were 
then  processed  using  the  SSOP  technique  as  in  Ref.  14. 
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Fig.  2  Comparison  of  normalized  images  taken  of  the  spatial  pattern  when  the  LEDs  were  individually 
illuminated  without  temporal  modulation  (top  row)  versus  temporally  demodulated  images  of  each 
wavelength  using  the  TM-SFDI  system  (bottom  row).  The  images  closely  match,  showing  the  ability  of 
temporal  demodulation  to  separate  different  wavelengths.  The  scale  bar  applies  to  all  images. 


obtained  via  Hilbert  transform  followed  by  inverse  Fourier 
transform  of  the  high  frequency  content.  The  demodulated 
AC  and  DC  images  were  compared  with  a  reference  standard 
to  obtain  the  calibrated  reflectance  of  each  pixel  in  the 
image.  Optical  properties  were  calculated  by  interpolating  a 
look-up  table  relating  calibrated  reflectance  values  to  optical 
properties,  which  was  generated  by  scaling  of  a  single  Monte 
Carlo  simulation.12,18 

2.3  Accuracy  and  Drift 

To  test  the  accuracy  of  the  TM-SFDI  system,  a  set  of  13  tissue 
mimicking  phantoms  with  varying  optical  properties  were 
imaged  three  times  over  3  days  with  a  commercial  SFDI  system 
(Oxlmager  RS,  Modulated  Imaging,  Inc.)  at  wavelengths  of 
526,  659,  and  729  nm  to  establish  the  “gold  standard”  measure¬ 
ments  of  the  phantom  optical  properties.  The  same  set  of  phan¬ 
toms  was  also  imaged  three  times  over  3  days  with  the  TM-SFDI 
system  at  wavelengths  of  519,  652,  and  740  nm.  Both  systems 
used  measurements  of  a  reference  standard  with  known  optical 
properties  for  instrument  response  calibration.  To  account  for 
variability  in  the  TM-SFDI  system,  the  reference  standard  was 
imaged  five  times  per  day  and  the  temporally  demodulated 
images  were  averaged  together  before  being  used  to  calculate 
//a  and  at  each  wavelength.  Optical  properties  of  the  phantoms 
from  both  the  commercial  system  and  the  TM-SFDI  system 
were  calculated  using  the  same  SSOP  algorithm  and  compared 
against  each  other.  Because  the  same  wavelengths  were  not 
available  for  both  systems,  a  correction  factor  was  imple¬ 
mented  to  estimate  the  optical  properties  calculated  by  the 
commercial  system  at  the  wavelengths  used  in  the  TM-SFDI 
system.  To  construct  the  correction  factor  for  //a,  the  absorption 
spectrum  of  nigrosin,  the  chromophore  used  to  make  the  tissue 
mimicking  phantoms,  was  measured  using  a  spectrophotometer. 


The  correction  factor  for  absorption  was  defined  as 
Aa (-^tm )  —  where  ^AATM)  is  the  estimated  absorp- 

tion  coefficient  at  the  TM-SFDI  wavelength  2XM,  Aa(^Com)  is  the 
measured  absorption  coefficient  at  the  wavelength  of  the 
commercial  system  (icomX  and  e(A)  is  the  extinction  coefficient 
of  nigrosin  at  the  wavelength  A  measured  with  the  spectropho¬ 
tometer.  To  estimate  how  //'  changes  with  wavelength, 
we  assumed  a  Rayleigh  scattering  (//s'  cx  A~ 4)  relationship  and 

defined  the  correction  factor  as  As(^tm)  =  ^s^cT4^Com,  where 

As(^Com)  is  the  reduced  scattering  coefficient  calculated  by 
the  commercial  system  at  ACom.  Rayleigh  scattering  was 
assumed  because  scattering  in  the  phantoms  is  controlled  by 
altering  the  concentration  of  titanium  dioxide  nanoparticles, 
the  diameters  of  which  are  much  smaller  than  the  wavelengths 
of  light  used.  The  maximum  change  in  optical  properties 
using  these  correction  factors  was  slightly  over  5%.  Each 
TM-SFDI  scan  consisted  of  128  2048  x  1536  pixel  images, 
with  an  FoV  of  8.5  X  6.4  cm.  To  improve  the  signal-to-noise 
ratio,  the  temporally  demodulated  images  were  downsampled 
using  2x2  averaging. 

To  assess  the  stability  of  the  TM-SFDI  system,  a  single  tissue 
mimicking  phantom  was  scanned  every  10  min  for  2  h  following 
a  30-min  warm  up  period,  and  the  variations  in  calibrated  reflec¬ 
tance  and  optical  properties  at  each  wavelength  over  that  time 
were  calculated  based  on  comparison  with  a  reference  standard 
imaged  at  the  start  and  end  of  the  drift  test.  Each  scan  consisted 
of  128  images. 

2.4  Effect  of  Ambient  Light 

To  test  the  effect  of  ambient  light  on  TM-SFDI,  a  calibration 
phantom  was  imaged  in  complete  darkness.  TM-SFDI  videos 
were  then  acquired  at  various  levels  of  ambient  light  controlled 
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by  changing  the  distance  between  the  bulb  of  a  standard  desk 
lamp  and  the  sample.  Images  of  a  single  wavelength  not 
modulated  in  time  were  also  acquired  to  compare  the  TM-SFDI 
technique  with  the  SSOP  method. 

3  Results  and  Discussion 

3.1  Accuracy 

We  found  that  there  was  a  good  agreement  in  optical  property 
extractions  between  the  gold  standard  commercial  instrument 
and  the  proof-of-concept  TM-SFDI  system  (Fig.  3).  The  average 
error  in  absorption  between  the  TM-SFDI  and  the  gold  standard 
was  13%.  The  TM-SFDI  system  was  more  accurate  at  extracting 
reduced  scattering,  with  an  average  error  of  8%  compared  with 
the  gold  standard.  The  absorption  measurements  were  biased, 
consistently  underestimating  the  actual  absorption  coefficient 
by  an  average  of  0.003  mm-1,  while  the  scattering  measure¬ 
ments  were  biased  in  the  other  direction,  overestimating  //' 
by  an  average  of  0.088  mm-1. 

In  addition  to  homogeneous  phantoms,  a  more  complex  tis¬ 
sue  simulating  phantom  was  imaged.  The  sample  was  a  3-D 
printed  block  of  acrylonitrile  butadiene  styrene  that  featured 
a  network  of  vascular  like  channels  filled  with  a  highly  absorb¬ 
ing  acrylic  copolymer  (P400SR,  Stratasis,  Inc.).  The  complex 
phantom  was  imaged  with  the  TM-SFDI  system,  and  a  good 
correlation  between  the  printed  channel  and  the  absorption 
map  was  noted  (Fig.  4). 

3.2  Drift 

The  drift  of  the  TM-SFDI  system  was  tested  by  allowing  the 
system  to  warm  up  for  30  min  and  then  repeatedly  measuring 
the  same  phantom  every  10  min  for  2  h.  The  average  changes  in 
calibrated  reflectance  over  this  time  period  were  less  than  3%  for 
the  AC  and  DC  components.  However,  these  errors  translated 
into  a  relatively  large  average  variation  in  pid  of  13%  about 
the  mean.  The  average  variation  in  yws'  was  1%.  There  was  no 


systematic  rise  or  fall  to  the  change  but  rather  a  random  variation 
about  the  mean.  To  determine  the  source  of  the  drift,  we  tested 
each  LED  individually  with  DC  illumination  and  no  temporal 
modulation  and  found  a  variation  of  1%  to  3%  of  output  inten¬ 
sity  over  time.  Though  relatively  small,  this  variation  would  also 
change  the  calibrated  reflectance  by  the  same  1%  to  3%,  which 
maps  to  a  greater  than  10%  change  in  fid.  The  observed  variation 
in  optical  properties  is  likely  due  to  instability  of  the  LEDs  and 
driver  electronics,  which  were  both  consumer  grade.  We  expect 
that  upgrading  these  components  and  implementing  temperature 
control  would  reduce  the  variability  in  measured  optical 
properties. 

3.3  Effect  of  Ambient  Light 

We  found  that  the  TM-SFDI  system  was  almost  entirely 
insensitive  to  the  ambient  light  level.  In  Fig.  5,  the  calibrated 
reflectance  for  planar  illumination  is  shown  for  SSOP  (blue) 
and  TM-SFDI  (red)  along  with  linear  fits  to  the  data  at  a  variety 
of  ambient  light  levels.  The  slope  of  the  linear  fit  of  the 
TM-SFDI  DC  reflectance  was  not  significantly  different  from 
zero  ( p  =  0.5)  using  a  two  sided  /-test,  indicating  that  there 
was  no  relationship  between  ambient  light  level  and  calibrated 
reflectance.  The  change  in  calibrated  reflectance  using  SSOP 
would  result  in  extremely  large  errors  in  optical  properties  as 
ambient  illumination  changed.  Future  instruments  using  TM- 
SFDI  will  need  to  avoid  frequencies  near  the  AC  power  supply 
(typically  50  or  60  Hz)  to  avoid  contamination  from  the  flick¬ 
ering  of  room  lights. 

3.4  Discussion 

The  TM-SFDI  system  presented  here  allows  for  the  rapid 
acquisition  of  diffuse  optical  images.  At  56  frames /s  and 
a  typical  video  length  of  128  frames,  we  were  able  to  capture 
spectroscopic  images  of  optical  properties  at  three  wavelengths 
in  2.3  s.  The  number  of  simultaneous  wavelengths  that  can  be 
collected  is  limited  by  the  number  of  video  frames  acquired.  The 
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Fig.  3  Comparison  of  TM-SFDI  with  gold  standard  in  terms  of  (a)  scattering  and  (b)  absorption.  The 
identity  line  is  shown  in  gray.  Error  bars  indicate  ±1  standard  deviation  of  the  average  measured  optical 
property  across  the  3  days.  Acquisition  time  was  2.3  s  for  the  TM  system  versus  ~30  s  for  the  gold 
standard  device. 
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Fig.  4  (a)  Photograph  of  the  3-D  printed  phantom  with  vessel-like  inclusions  under  strong  back  illumi¬ 
nation.  (b)  Temporally  and  spatially  demodulated  image  of  calibrated  DC  reflectance  at  652  nm.  (c)  Map 
of  the  absorption  coefficient  and  (d)  map  of  the  reduced  scattering  coefficient  of  the  sample  (d)  at  652  nm. 
Scale  bar  applies  to  all  images.  Median  filtering  was  used  on  panels  (c)  and  (d)  to  remove  artifacts. 


Insensitivity  to  background  light 


Fig.  5  Temporally  modulated  (blue)  SFDI  compared  with  single  shot 
SFDI  (red)  under  varying  levels  of  ambient  light.  TM-SFDI  is  largely 
insensitive  to  the  ambient  light  level  while  SSOP  is  highly  sensitive. 
Error  bars  represent  ±1  standard  deviation  of  the  pixels  in  the  image. 
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maximum  number  of  wavelengths  that  can  be  unambiguously 
demodulated  is  no  more  than  one  per  FFT  bin  (excluding 
DC).  This  means  that  a  128  frame  scan  could  theoretically  sup¬ 
port  63  wavelengths,  while  a  256  frame  scan  could  support  127, 
although  practically,  this  number  of  wavelengths  may  be 
reduced  due  to  spectral  leakage.  We  experimented  with  using 
fewer  video  frames  (64,  32,  16,  etc.)  and  found  an  increased 
scan-to-scan  variability.  Using  64  images,  there  were  only  mod¬ 
estly  higher  errors  (10%  in  scattering  and  24%  in  absorption), 
meaning  acquisition  speeds  of  around  1  s  with  reasonable  errors 
are  possible.  Since  the  maximum  frame  rate  of  the  camera 
increases  when  fewer  CCD  lines  need  to  be  read,  using  a  smaller 
FoV  could  also  increase  acquisition  speed.  For  an  image  of 
512x512  pixels,  it  would  be  possible  to  increase  the  frame 
rate  from  55.6  to  80  fps,  resulting  in  a  minimum  acquisition 
time  of  800  ms  for  up  to  16  wavelengths. 

Recent  work  has  utilized  a  similar  temporal  modulation  of 
excitation  wavelengths  for  single  pixel  SFDI.19  Despite  advan¬ 
tages  in  simplicity  of  hardware,  cost,  and  bandwidth,  single 
pixel  imaging  scales  poorly  with  resolution,  requiring  the  pro¬ 
jection  of  over  31,000  patterns  for  each  spatial  frequency  and 
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phase  to  acquire  an  image  with  a  comparable  number  of  pixels 
used  in  this  study.  For  example,  Torabzadeh  et  al.19  used  single¬ 
pixel  SFDI  to  acquire  400  samples  in  80  ms,  meaning  a  single 
high  resolution  image  could  be  collected  every  6.2  s.  While  a 
large  number  of  pixels  is  not  strictly  necessary  for  low-resolu¬ 
tion  diffuse  imaging,  it  allows  for  pixel  binning  to  reduce 
the  overall  noise  level.  Greater  numbers  of  pixels  would  also 
allow  for  the  possibility  of  imaging  larger  FoVs  without 
impacting  resolution. 

The  TM-SFDI  system  has  several  limitations  worth  noting. 
The  first  is  variability  in  extractions.  The  errors  in  absorption 
are  likely  due  to  the  fact  that,  at  the  chosen  spatial  frequencies, 
small  changes  in  reflectance  map  to  large  changes  in  pa,  mean¬ 
ing  that  the  effects  of  noise  and  LED  instability  are  amplified.4 
It  is  likely  that  temperature  variations  of  the  LEDs  and  driver 
circuit  led  to  instability  of  the  LED  output  and  resulted  in 
uncertainty  in  the  absorption  estimates.  We  also  found  that 
the  estimates  for  optical  properties  were  biased  with  the 
absorption  being  underestimated  and  the  scattering  overesti¬ 
mated.  We  believe  both  biases  are  due  to  broadband  noise  con¬ 
taminating  the  images,  resulting  in  small,  systematic  errors  in 
diffuse  reflectance  as  well  as  nonlinearities  in  system  response, 
which  are  not  accounted  for  in  the  calibration  step  and  result  in 
increased  errors  for  samples  with  optical  properties  that  differ 
significantly  from  the  calibration  phantom.  It  is  likely  that 
improvements  in  the  signal-to-noise  ratio  and  careful  choice 
of  calibration  phantom  will  help  resolve  this  issue.  Second, 
the  acquisition  time  of  1  to  2  is  relatively  long  compared 
with  the  SSOP  method,  and  the  sample  must  be  held  motionless 
for  that  entire  time  to  achieve  a  clear  image  of  the  resulting 
optical  properties.  This  could  potentially  be  solved  by  adding 
an  image  registration  step  to  align  the  frames  of  the  video 
prior  to  processing.  Third,  although  scans  can  be  performed 
relatively  quickly,  processing  of  the  data  requires  significant 
computational  resources.  Each  scan  requires  over  300,000  FFTs 
to  be  performed  as  well  as  searching  a  nonlinear  look-up-table 
for  each  pixel  at  each  wavelength.  On  a  relatively  modern  desk¬ 
top  computer  (Intel  i7-6700,  16  GB  RAM),  the  time  to  process 
each  scan  was  about  80  s.  However,  much  of  the  processing  is 
parallelizable  and  could  likely  be  vastly  sped  up  by  offloading  to 
a  graphics  processing  unit. 

4  Conclusion 

In  conclusion,  we  have  constructed  a  high-speed  SFDI  system 
that  enables  the  acquisition  of  spectrally  resolved  optical  proper¬ 
ties  without  needing  to  switch  between  illumination  patterns 
or  wavelengths.  Currently,  the  system  operates  at  ~750  ms/ 
wavelength,  but,  since  acquisition  time  does  not  scale  with 
the  number  of  wavelengths  used,  this  could  be  increased  to 
~50  ms/ wavelength  by  adding  additional  illumination  wave¬ 
lengths.  By  mapping  the  spectral  components  of  light  into  the 
time  domain,  this  system  eliminates  the  need  to  acquire  each 
wavelength  sequentially.  Importantly,  the  temporal  modulation 
of  light  also  makes  this  system  insensitive  to  background  illu¬ 
mination,  making  it  a  practical  technique  for  translation  into  the 
clinic,  where  light  levels  vary  widely  and  cannot  easily  be  con¬ 
trolled.  Future  work  will  investigate  increasing  the  stability  of 
the  LEDs,  increasing  the  number  of  modulated  wavelengths 
using  dispersed  supercontinuum  light  in  conjunction  with  a  lin¬ 
ear  spatial  light  modulator,  testing  the  TM-SFDI  on  a  variety  of 
biological  samples,  and  including  single-shot  height  correction 
for  more  accurate  imaging  of  3-D  samples.15 
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